Mechanically induced metal— insulator transition in carbyne 
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First-principles calculations for carbyne under strain predict that the Peierls transition from 
symmetric cumulene to broken-symmetry polyyne structure is enhanced as the material is stretched. 
Interpretation within a simple and instructive analytical model suggests that this behavior is valid 
for arbitrary ID metals. Further, numerical calculations of the anharmonic quantum vibrational 
structure of carbyne show that zero-point atomic vibrations alone eliminate the Peierls distortion in 
a mechanically free chain, preserving the cumulene symmetry. The emergence and increase of Peierls 
dimerization under tension then implies a qualitative transition between the two forms, which our 
computations place around 3% strain, whereupon the conductivity of carbyne sharply decreases. 
Thus, zero-point vibrations and mechanical strain jointly produce a change in symmetry resulting 
in the transition from metallic to insulating state. In any practical realization, it is important 
that the effect is also chemically modulated by the appropriate choice of terminating groups. Our 
findings are promising for applications such as electromechanical switching and band gap tuning via 
strain, and besides carbyne itself, they directly extend to numerous other systems that show Peierls 
distortion. 



Carbyne — the linear allotrope of carbon — is perhaps 
one of the most unusual materials due to its ultimate one- 
atom thinness. Although this elusive material is notori- 
ously hard to prepare and has been perceived as exotic or 
even completely fictitious, the development of methods to 
synthesize it proceeds at a steady rate, with input from 
both experimentalists and theoreticians [1-7J. Among 
the most notable recent achievements, chains with length 
of up to 44 atoms [8j and such complex molecular mech- 
anisms as carbyne-based rotaxanes have been syn- 
thesized. This progress is driven by carbyne's attractive 
physical properties such as unusual electrical [11] [12] and 
mechanical behavior [13j, or its large specific area [l4j. 
Accordingly, a better theoretical understanding of this 
material is becoming more and more relevant. 

It has long ago been established by the quantum chem- 
istry community that carbyne undergoes the Peierls 
transition p!6l4T9] that converts it from the cumulene 
(=C=C=)n to the polyyne (-C=C-)n form. More re- 
cently, it has been suggested that the zero-point (ZP) 
vibrations may be able to overcome the Peierls instabil- 
ity and restore the nonalternating bond pattern in car- 
byne [20j, as well as in another canonical example of 
Peierls instability, polyacetylene [21j. Since the symmet- 
ric and broken- symmetry forms have very distinct elec- 
tronic properties (metallic and semiconducting, respec- 
tively), investigations of this effect are important from 
not just the fundamental chemical physics perspective, 
but also practically relevant for the field of conducting 
polymers. 
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Recently, we have found through first-principles cal- 
culations that stretching carbyne produces a number of 
interesting consequences [13j (also observed experimen- 
tally after the present study had been completed [22j). 
The bond length alternation (BLA), defined as the dif- 
ference between the single and triple bond length, the en- 
ergy of Peierls distortion, and the band gap all increase 
with strain. The latter is counterintuitive because nor- 
mally one would expect strain to result in a weaker or- 
bital overlap, and hence, smaller tt — tt* band splitting. 
The BLA increase could counter this effect via symmetry 
breaking-induced increase of this splitting. However, the 
BLA increase itself can not be easily explained within a 
classical 'balls-and-springs' mechanical representation of 
carbyne. 

Thus, we begin with an exploration of the effect of 
strain on the Peierls transition. We use a simple analyt- 
ical model to explain why strain amplifies the Peierls ef- 
fect. Then, we investigate through first-principles Born- 
Oppenheimer potential energy surface calculations how 
the atomic vibrations manifest themselves in equilibrium 
and under strain, confirming the earlier predictions that 
the Peierls effect at zero strain is too weak, so that the 
transition barrier lies well below the ZP vibrational level. 
Since this barrier is increased by strain, we arrive at the 
prediction of a strain threshold whereupon carbyne tran- 
sitions from the ZP vibrations- to Peierls instability- 
dominated regime, accompanied by a sharp change in 
electrical conductivity. Our first-principles numerical 
analysis validates this prediction, placing this transition 
around ~3% strain (1.5-5.5%). Besides carbyne, our 
findings naturally suggest that this novel physical effect 
must also be observable in other one-dimensional systems 
showing Peierls-transition behavior, such as conducting 
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FIG. 1: The (a) bond length alternation and (b) band gap 
of carbyne as a function of strain using the PBE and HSE06 
density functionals. The inset in (a) shows the polyyne- 
cumulene energy difference as a function of strain (Peierls 
barrier) . 
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polymers or carbon nanotubes [23j. 

We begin by presenting our first-principles data show- 
ing the amplification of the Peierls instability under ten- 
sion. Fig. [l] (a) plots the BLA as a function of strain, 
and the corresponding energy difference between alter- 
nating (polyyne) and nonalternating (cumulene) struc- 
tures is plotted in the inset. The band gap dependence 
on strain is shown in Fig. [l](b). The Peierls distortion 
is known to be extremely sensitive to the electronic ex- 
change interaction and, in particular, poorly described by 
regular density functional theory [24j, therefore the cal- 
culations were repeated using two methods: (blue) the 
PBE generalized-gradient density functional [25l |26] and 
(purple) the HSE06 hybrid functional that includes exact 
exchange [27^ 28j. All Fig. [l] plots conspicuously demon- 
strate the same universal behavior of increase with strain 
(also observed in high-level GW calculations [13l [22]). 

Clearly, the effect observed so pronouncedly in Fig. [l] 
data calls for an explanation. We can understand the 
behavior of Peierls transition under strain with the aid 
of the following simple Kronig-Penney type model. Con- 
sider a chain with two atoms per period a with coor- 
dinates xi = and X2 = (a — 6) /2, where b is the 
BLA. The atoms are represented as attractive Dirac 
delta functions, and the potential in the unit cell is 
V (x) = —S{x — Xi) — S{x — X2). The textbook expres- 
sion for the electronic part of the Peierls effect is [29j 
Ee (X log \ Vk \ /hk^ where k is the Brillouin zone edge 
{k = ir/a) and Vk is the corresponding Fourier component 
of the lattice potential. For our delta-function model, it 
is simply Vk = sin {nb/a). 

The above electronic contribution is counteracted by 
the lattice deformation energy, Ei = Cb'^ + O (b^) . The 
spring constant C generally depends on a in some un- 
known way. Clearly, C should decrease with strain. 
For mathematical simplicity, we choose C = Co /a. 
Now all that remains is to minimize the total energy 
with respect to the BLA 6, the condition for which is 
d {Eq El) jdb = 0. But, upon examination, one finds 
that with our choice of model, the derivatives of both 



FIG. 2: Analytical model for Peierls instability using a = 
(yellow) , (green) 1.0, (purple) 1.1 of the equilibrium 
lattice constant ao [arbitrary units] and two models for the 
lattice deformation potential: strain-dependent-stiffness har- 
monic spring as in the text ("H", dotted lines), constant- 
stiffness harmonic spring (black line) , and the Morse potential 
(dashed lines, "M"). Panel (a) shows the graphical solution for 
the balance of Peierls (Fe) and restoring forces, and the cor- 
responding Peierls instability energies are shown in (b). 



energy terms are only functions of the bja ratio. This 
immediately means that b/a = const ^ and hence, b oc a; 
that is, the BLA increases with strain. 

To demonstrate the robustness of our conclusion with 
respect to changes in the C (a) dependence. Fig. [2] (a) 
illustrates the graphical solution for the balance of the 
electronic driving force of Peierls instability (solid curves 
colored according to the strain) and the restoring forces, 
using the harmonic C = Co/a model (dotted lines, "H"), 
a constant spring stiffness Co /do (solid black line), and 
the Morse potential (dashed lines, "M"). Evidently, the 
effect of BLA increase with strain should hold for essen- 
tially any reasonable form of the restoring force poten- 
tial. Finally, Fig. [5] (b) plots the corresponding energy 
differences between polyyne and cumulene forms, demon- 
strating that our model can indeed explain all three plots 
of Fig. [l] (We note that the Morse potential, due to its 
rapid softening with length, can produce a nonmonotonic 
dependence of distortion energy at large tensile strain — 
see dashed line in Fig. [2] (b).) 

Up until now, we have confined ourselves to the static 
picture of carbon atoms. Within this picture, only the 
cumulene form of carbyne is stable (at low temperatures) , 
and stretching it produces merely quantitative changes in 
the band gap — a phenomenon known previously in solid 
state physics. However, as we demonstrate below, the 
inclusion of quantized nuclear motion changes the pic- 
ture dramatically. If ZP vibrations can overpower the 
Peierls instability in mechanical equilibrium and restore 
the symmetric metallic structure as previously suggested 
[2Q[ [21], then, considering this together with our above 
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findings, one is led to expect ttiat at some critical strain 
value, the stronger Peierls instability may reverse this 
phenomenon — causing a qualitative change in properties. 
Our calculations indeed confirm this. 

For the vibrational analysis, the Born-Oppenheimer 
potential energy surface scans along the BLA coordinate 
were computed with the PBE and HSE06 functionals, 
yielding the potential energy surface (PES) profiles Ea (b) 
as the lattice constant a was increased from its equilib- 
rium value. Because of the strong anharmonicity of the 
"inverse McDonald's^^" cj-shaped PES profiles, we had 
to turn to numerical calculations of the vibrational levels 
at each strain value. We utilized the Fourier grid Hamil- 
tonian method [3Ql [3T] which had been previously suc- 
cessfully used to study polyacetylene in a similar context 
[21j. Following the previous studies, we assume that the 
longitudinal vibrations can be decoupled from the two 
transverse phonons and treated separately, which is jus- 
tified by the non-negligible bending stiffness of carbyne 
[13l [32H34] . (Indeed, the inclusion of transverse vibra- 
tions and dispersion does not change our conclusions, as 
detailed in the Supplemetary Materials.) We used a uni- 
form grid of 250 points in —0.4 A < b < 0.4 A, using in- 
terpolation of DFT-computed values of energy, and 6.005 
a.m.u. — the reduced mass of two carbon atoms — for the 
oscillator mass. 

The results are presented in Fig. [3j which is the central 
figure of this paper. The red lines show the cj-shaped PES 
(in eV), blue lines denote the computed vibrational lev- 
els with filling below the zeroth level, and the green lines 
plot the ZP wavefunction, as computed using the HSE06 
functional. The right panels show the reconstructed real- 
space densities of nuclei at the respective strain values. 
The profile along the horizontal axis is taken directly 
from the vibrational wavefunctions, and an arbitrarily 
chosen value of 0.12 A is used for the transverse broad- 
ening. 

Clearly, the mechanically relaxed system has a single 
wavefunction maximum, corresponding to the cumulene 
structure of carbyne — although the PES has a local maxi- 
mum in this point. In accordance with the above discus- 
sion, the maximum ('Peierls barrier') rises with strain, 
and at 3%, the wavefunction already has a slight dip in 
the center, resulting in a picture of 'elongated' atoms. At 
7% strain, the ground-state wavefunction is clearly sep- 
arated in two blobs, and the first excited level (antisym- 
metric with a node in the middle) approaches the ground- 
state level, indicative of the transition to the double-well 
potential-like regime. Finally, at 10% strain, the first 
excited level is well below room-temperature ksT rela- 
tive to the ground state, as in two independent potential 
wells. The second level dives below the Peierls barrier, 
and the third begins to approach it. (At 15% strain these 
two also become degenerate; data not shown.) 

The above-described sequence of events has a very pro- 
found consequence. As the strain is increased, carbyne 
transitions from the conducting electronic behavior of cu- 
mulene to the nonconductive polyyne structure with its 
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FIG. 3: (left) Evolution of the vibrational level structure of 
carbyne with strain: from 0% to 10%. Red lines represent 
the cj-shaped potential energy as a function of BLA, calcu- 
lated with the HSE06 density functional. Blue lines show the 
vibrational levels (the filled region is below the ZP level). 
Green lines show the ZP wavefunction. (right) Real-space 
atomic density distributions calculated based on the ZP wave- 
function. The plot in the bottom shows the strain dependence 
of Peierls barrier (circles, PBE; squares, HSE06) and vibra- 
tional levels (colored x and diamond symbols). 



sizable band gap (see Fig. [T] (b)). The exact point at 
which this happens is not clear from Fig. [3] plots, and di- 
rect phonon-coupled electronic transport calculations are 
infeasible. However, an approximate criterion can be sug- 
gested basing on the bottom plot in Fig. [3) which shows 
the joint evolution of the Peierls barrier and vibrational 
levels under strain. The wavefunction starts to acquire 
the dual-maximum character when the ground-state level 
dives below the Peierls barrier, and we can consider the 
transition to be well over by the time the same happens 
to the first excited level. As seen in the bottom of Fig. [3| 
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FIG. 4: Effect of endgroups on the BLA in carbyne. Termi- 
nation of the chain with sp^ groups (methyl — CH3 and phenyl 
— CeHs) increases the BLA, whereas the sp^ methylene group 
(— CH2) decreases it, as compared to an unterminated infinite 
chain (black dashed line). The data for 12- and .uu-t^^^ 
chains show weak decrease of BLA with chain length away 
from the ends with methylene termination. 

these two events happen at the 1.5% and 5.5% strain val- 
ues, respectively, which is a reasonably mild strain from 
the mechanical point of view. (PEE gives 8% and slightly 
above 10%, correspondingly.) 

So far our calculations treated carbyne as an infinite 
chain. In a real experiment, however, the ends of finite 
carbyne chains inevitably represent special points that 
might affect the picture laid out above. (Note that our 
predictions should still apply directly to carbyne rings.) 
The PES of infinite carbyne contains pairs of minima and 
is invariant under translation by half a unit cell (assum- 
ing polyyne form), which corresponds to the exchange of 
triple and single bonds. But, as a simple example, the 
two finite structures described by chemical formulae 

H - C = C (-C = C-)^ C = C - H 

and 

H = C- C(=C-C=)^C-C = H 

are, clearly, very nondegenerate, because hydrogen can- 
not form stable triple bonds with carbon. On the other 
hand, the use of end-groups that locally favor the cumu- 
lene structure, such as methylene =CIl2, could be ex- 
pected to be less disruptive. The effect of end-groups 
will clearly have a different strength for different moi- 
eties, and a comprehensive investigation of the variety of 
possibilities is an extensive task (further aggravated by 
complexity of collective atomic motions in such many- 
atom systems) that is beyond the scope of the present 
study. However, simply looking at how the BLA depends 
on the endgroup type in finite carbyne fragments apready 
provides a useful glance into the issue. 



The results of our finite-chain calculations with dif- 
ferent terminating groups are presented in Fig. |4] The 
strip from 1.27 to 1.30 A denotes the bond lengths in the 
ideal infinite polyyne chain (ignoring the ZPV). The bond 
lengths in the 5p^-terminated systems (-CH3, -CqR^) 
fall outside this strip, showing the enhanced polyyne-like 
character of the system. On the other hand, 5p^-type 
termination with -CII2 reduces the BLA effectively to 
zero (< 0.01 A for a 12-atom chain and < 0.005 A in 
the middle of a 100-atom chain), suggesting that this 
type of passivation does not significantly disrupt the po- 
tential energy picture shown in Fig. [3j Thus, the use 
of end-groups that locally favor the cumulene structure 
should allow our predictions to be observed. (Indeed, 
such groups are preferrable from the electronics applica- 
tions standpoint in order to provide metallic 'leads' at 
the ends of the chain.) In summary, we see that in ex- 
perimental realizations, the effect may be destroyed by 
improper chemical termination of carbyne chains, and 
thus, precise control over both the tension and chemistry 
is essential. 

Our findings regarding the importance of vibrations 
in carbyne add another example of a system where quan- 
tum effects in the motion of atoms have big consequences. 
This notably contrasts with the physicist's intuition say- 
ing that the carbon atom is too heavy for quantum effects 
to be significant, yet its mass is compensated for by the 
extreme mechanical stiffness of the carbyne lattice |13| . 
Carbyne represents one more case where vibrations, even 
at zero temperature, may stabilize a material in a con- 
figuration that is not even a local minimum of the Born- 
Oppenheimer PES [35]. Recently, quantum- mechanical 
effects in the nuclear motion were demonstrated to play 
roles in the dynamics of such systems as vacancies in 
graphene [36j or dislocations in 3D crystals [37J. How- 
ever, to our knowledge, carbyne (and likely, by analogy, 
polyacetylene) is the first example where such effects can 
lead to qualitative differences in the material's behavior — 
here, electronic conductivity. 

In conclusion, we have explored how the Peierls in- 
stability of linear carbon chains is enhanced by strain 
using first-principles calculations, and explained it with 
a simple analytical model. At zero strain, the Peierls 
instability in carbyne is too weak to overcome the quan- 
tum zero-point vibrations, and the stable structure cor- 
responds to the metallic cumulene rather than insulating 
polyyne — despite the lower potential energy of the latter 
structure. However, the enhancement of Peierls instabil- 
ity results in the reversal of this effect, and carbyne starts 
to demonstrate the polyyne-like character of the nuclear 
wavefunction, switching from the conducting to insulat- 
ing behavior around ~3% strain. This novel effect natu- 
rally applies to carbyne rings and other one-dimensional 
systems such as polyacetylene or carbon nanotubes. Our 
findings are important both as an interesting fundamen- 
tal physical effect and highly practically relevant for the 
science of conducting polymers. Although at present car- 
byne remains a highly exotic material, its unique proper- 
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ties will continue to fuel the community effort to achieve 
its synthesis in practically useful quantities. 

The computations were performed using the VASP 
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